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Abstract 

We generalize the fractal bubble model (FB), recently proposed in 
the literature as an alternative to the standard ACDM cosmology, to 
include a non-zero cosmological constant. We retain the same vol- 
ume partition of voids and walls as the original FB model, and the 
same matching conditions for null geodesies, but do not include ef- 
fects associated with a nonuniform time flow arising from differences 
of quasilocal gravitational energy that may arise in the coarse-graining 
process. The Buchert equations are written and partially integrated 
and the asymptotic behaviour of the solutions is given. For a universe 
with A = 0, as it is the case in the FB model, an initial void fraction 
with hyperbolic curvature evolves in such a way that it asymptotically 
fills completely our particle horizon. Conversely, in presence of a non 
vanishing A, we show that this does not happen and the voids fill a 
finite fraction f Voo < 1, where the value of (1 — / Uoo ) is expected to 
depend on A and the initial fraction f vi and also to be small. For its 
determination, a numerical integration of the equations is necessary. 
Finally, an interesting prediction of our model is a formula giving a 
minimum allowed value of present day dark energy as a function of the 
age of the universe and of the matter and curvature density parameters 
at our time. 

PACS numbers: 98.80.-k, 95.36.+X, 98.80.Es, 98.80.Jk 
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1 Introduction 



All the observations of the past decade (see in particular [U El [3j) are in 
agreement with an homogeneous and isotropic universe on large scales which 
seems to be currently in accelerating expansion. In the standard picture 
based on the Friedmann-Lemaitre models (FRLW), an accelerating universe 
invokes the presence of the so called dark energy. In the FRLW picture, 
this dark energy is given in terms of the cosmological constant A, leading 
to the standard cosmological model, i.e. the ACDM one. As a result, about 
70% of the present universe would be composed of an undetectable non- 
local vacuum energy. This dark energy represents a puzzle and perhaps 
the biggest problem in modern cosmology. In fact, a direct detection of a 
cosmological constant is still lacking. However, the dark energy is only the 
outcome of the Friedmann assumptions (homogeneity and isotropy) justified 
by the Copernican principle. In this respect, the observed inhomogeneities at 
scales < 200 Mpc give us the possibility to explore the role of inhomogeneity 
in cosmological models. In the last decade many attempts have been made 
(see for example [4]- [27]) to build physically sensible inhomogeneous models. 

One proposed approach (see for example [31 El El El El Uni HT] ) showed that 
large inhomogeneities can mimic an accelerating universe by using Lemaitre- 
Tolman-Bondi (LTB) spacetimes [281 EH], although only under successful 
completion of several conditions [12] . In this case, the cosmological constant 
can be mimicked only by supposing that a huge spherically symmetric giant 
void (or hump, see [30]) of the order of some Gpc is present, a rather unlikely 
assumption which is not in agreement with the CMB data (see also |31|). 

A different approach is averaging the spatial inhomogeneities (see |13j - 
|23j). This approach fulfils the Copernican principle, although in a statistical 
sense. Hence, the Copernican cosmologists introduce inhomogeneities with- 
out symmetries and then try to understand the modifications to the average 
evolution from backreaction. So far, the backreaction from small scales could 
be not enough to rule out dark energy [32l [33] , but the debate on this issue 
is still open. In this context, Wiltshire [221 [23] recently proposed an inter- 
esting model, the 'fractal bubble' (FB) or 'timescape' cosmology [31], based 
on the Buchert average scheme [17] with two-scales mimicking the voids and 
walls in the observed web structure of our present day universe. The expan- 
sion history, usually attributed to a cosmological constant, is explained by 
a new physical effect: the different rate of clocks located in typical galaxies 
which are not expanding where the metric is approximately spatially flat, 
as compared to the clocks in expanding voids where spatial curvature is as- 
sumed to be negative. A uniform Hubble flow is imposed to solve the so 
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called Sandage-de Vaucouleurs paradox. It is the introduction of this gauge 
that allows clock effects. According with the estimates of |221 12,'jj . clocks of 
isotropic observers at a volume-average location in voids presently run at 
38% faster than the ones in walls. While this result is perhaps question- 
able, it is justified by Wiltshire as being the cumulative effect of a relative 
regional volume deceleration which is calculated [35] to be typically of order 
10~ 10 ms~ 2 for most of the life of the universe, well within the expectations 
of the weak field regime. 

In any case, although the cosmological constant remains a mystery for 
fundamental physics, the standard ACDM model remains the most convinc- 
ing cosmological model, with much astrophysical evidence in its support. 
The standard cosmological model relies on a universe homogeneous at any 
scale for any time t, a rather unrealistic assumption. 

As a consequence of the reasonings above, it is interesting to study the ef- 
fects of the observed inhomogeneous structure of the universe in presence 
of A, as for example the expansion of voids. The standard picture is to use 
perturbation theory to the standard ACDM model. To build an inhomoge- 
neous cosmological model with A and without perturbation theory we use 
the Buchert formalism. Within this formalism, we can gain advantage from 
the two-scale model obtained by Wiltshire, composed of walls and voids. 
However, remember that the model [221 [23] has been built to explain dark 
energy without A. The observed broadly uniform Hubble flow justified in 
[22l [23] the use of the uniform Hubble gauge. It is this gauge that allows a 
non- uniform time flow. In [22, 23j, the non- uniform time flow is described by 
the phenomenological lapse function 7. However, as we show in this paper, 
we can build our inhomogeneous model with A 7^ by taking the assump- 
tion of a uniform time flow (rather than a uniform Hubble flow) simply by 
changing in |22[l23j the interpretation of 7, which in our context is no longer 
a phenomenological lapse function. It is important to stress again that the 
clock effects depicted in [22[ [23] are only the outcome of the uniform Hubble 
flow used. However, note that no cosmological evidence of a non-uniform 
time flow has been observed. 

In this paper we use the mathematical structure of [22], [23] , but with A / 
in the Buchert equations and with a phenomenological lapse function set to 
unity. In this way, we have at our disposal a non perturbative model that 
will permit us to estimate, for example, the fraction of dark energy due to 
the observed inhomogeneities of our visible universe. Finally, we can study 
the evolution of voids in presence of A and compare this with future cos- 
mological data and with the same evolution obtained with models without 
A. 
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In section 2 we present the model by writing, step by step, the Buchert 
equations in a workable form. In section 3 we analyze the solutions of our 
model. Section 4 collects some final remarks and conclusions. Appendix is 
devoted to the study of the distance-redshift relation and to the matching 
conditions underlying our model. 



In this section we generalize step by step the calculations made in \22\ [23] by 
introducing a non-vanishing cosmological constant together with a uniform 
time flow gauge rather than the Wiltshire uniform Hubble flow gauge. As a 
result, clock effects disappear. Our aim is to build an inhomogeneous model 
within the cosmological constant paradigm which includes a description of 
the observed web structure of our visible universe and preserves the universal 
time flow of the concordance ACDM model. We want to mimic the observed 
web structure of the universe by means of a two-scale model. 

2.1 Buchert equations for the two-scale model 

As a first assumption, we suppose the existence of an almost uniform time 
flow. Hence, all the observers in our whole particle horizon, independently 
on their spatial location, use for time measurement the same coordinate 
t. In the Buchert formalism, we can introduce Gaussian coordinates (t,x l ) 
which are comoving with the fluid that is assumed to be dust filled and 
irrotational. The volume average is given to be our present particle horizon. 
For any scalar quantity il){t,x % ), the average with respect to the volume 
horizon % is: 



2 The model 



< 4>{t,x l ) > 



h — 



1 




X 




where denotes the determinant of the three metric on the slice a t 
const, and 





The dimensionless effective scale factor a%{£) is given by 




(3) 
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The averaged expansion rate < 9 >u is defined by: 

<0> H = ^ = ^, (4) 

where dot indicates the time derivative. For the Hubble flow we have H = 
< 2™ ■ In what follows we drop the subscript %. The relevant exact Buchert 
equations (see [17]) are (we use units with c = 1): 

3^- = 8nG < p > + A-l<K>, (5) 
a z 2 2 

3- = -4vrG < p > +Q + A, (6) 
a 

< p > = -3- < p >, (7) 
a 

Q = \[<d 2 >-<d> 2 } -2<a> 2 , (8) 

where Q is the kinematical backreaction, a represents the shear and 1Z is the 
curvature of the hypersurface t = constant. For the partitioning we have in 
mind, which is similar to the FB one, if we average over sufficiently large 
regions (see [35]), then the averaged shear < a > can be neglected. The 
integrability condition for the system (E])-([H]) is given by 

{<oQ + 2<1l>)a + a Q+ <TZ> = 0. (9) 

Equation ([5]), in terms of the density parameters, can be put in the form 
(the so called cosmic quartet) 

n m + n A + n k + n Q = 1, (10) 
n m = 87r 3 ^? > , n k = = -^p, n A = (ii) 

where H is the Hubble rate function. 

For the partitioning of our visible universe, we follow the two-scales 
model of the original FB model (see \22\ [23] for more). First of all we 
must specify the so called finite infinity regions. In the universe, there exist 
regions expanding and contracting. However, there are regions, called finite 
infinity Fj (see [23] [36] ) , that are the set of timelike boundaries of compact 
disjoint domains I for which the average expansion is zero and becomes 
positive outside, i.e. for any surface t = const: 

< 8 > UjFj = 0, (12) 
3 Dj, Fi D Di -> < 9 > U/D/ > 0. 
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In other words, {JjFj contains regions which on average are non-expanding. 
Our position is within a Fj region. The Fj regions are well within the wall 
regions, which are defined as the disjoint regions that have on average a 
spatially flat metric with expansion factor a w (t): 

ds 2 w = -dt 2 + a 2 w [dr& + vldfl 2 ] . (13) 

Outside the wall regions (the region complementary to the walls with respect 
to the particle horizon) we have the void ones that are defined as expanding 
regions where the metric is on average hyperbolic with scale factor a v (t) 

ds 2 v = -dt 2 + a 2 v [drj 2 + sinh 2 ^)^ 2 ] . (14) 

Another crucial assumption we make (also in the FB model) is the presence 
of a scale of homogeneity. This assumption is justified by the fact that 
statistical analysis of the galaxy distributions from available data seem to 
indicate an homogenization on scales of order of (100 — 200)Mpc. By scale 
of homogeneity we mean that averages of any variable beyond this scale 
practically will no longer depend on the scale. For the Hubble function in 
walls and voids we have: 

TT < > w a w < 8 > v d v . . 

H w = = — , H v = — - — = — . (15) 

With 'a' we denote the scale factor at the scale of homogeneity. The Buchert 
equations for a(t) are given by ©-([H]). Moreover, in what follows, we denote 
with the subscript T an initial early time ij. Since the main cosmological 
data (distance-redshift relation) are available starting from the recombina- 
tion era, we intend with ij this early time. Furthermore, we denote with 
Vi the initial volume of the particle horizon at the recombination era. For 
the average on the volume representing the whole visible universe (particle 
horizon) defined by V = a 3 Vi we have 

O — fwi Q"u) ~\~ fvi Q"ui fwi ~\~ fvi — 1; (^) 

where f w i and f„ are initial fractions of walls and voids. Formula (|16p relates 
the expansion factors at the two-scales (walls and voids) to the average 
geometry represented by a(t) at the scale of homogeneity. The wall volume 
is V w = V w id w while the void volume is V v = V v ia%, being f w i = and 
fvi = ^r- For a time t > t{ the picture is 

^3 ^3 
fv(t) ~\~ fw(t) = 1; fw = fwi 3") fv = fvi q- (1^) 



6 



At the volume horizon, for the Hubble rate we have 

H = f w H w + f v H v . (18) 
Moreover, we define I w and I v according to 

H = I W H W = I V H V . (19) 

In the Wiltshire model we have, instead of I w , the phenomenological lapse 
function 7. Differently to the FB model, in our model I W ,I V are not related 
to possible time delay effects. In practice, we adopt the same partitioning 
present in [23], but in a different context, i.e. with a uniform time flow gauge 
instead of a uniform Hubble one. As a consequence, in our context I w is 
only a measure of the ratio between the averaged Hubble flow at the scale 
of homogeneity and the one at the wall scale. This choice can be justified 
by the fact that, at present day, no evidence for a non-uniform time flow 
at cosmological scales has been observed, although the converse cannot be 
excluded. By derivating f v (t) given by (fTT|) we have: 

t. - -l - 3d - /,)d - o g - 3, i;(;%" )H . (20) 

where h is defined by I v = hl w . For a general scalar -ip we have: 

<1p >= fwi>w + fvlpv (21) 

As a consequence of (f2Tj) we get: 

< e 2 >= 9(f w Hl + f v H 2 v ). (22) 



From (j8]). (118p . (I20p and (J22J) and neglecting the possibility of intrinsic back- 
reaction within the walls and voids (see|20} [37]) we obtain: 

n af r-i t\iu u \2 6/y(l — /d)(1 - h) 2 H 

Q = 6/„(l - fv){H v - H w ) = — — — ——2 — • (23) 

[h + (l- h)f v \ 



A useful alternative expression for the back-reaction Q is: 



2f 2 

Q= — • (24) 

^ 3/41-/,) 1 > 



Concerning the mean three curvature < 1Z >, we have < 1Z > v = 6% and 
< 1Z > w = 6% = since by construction k w = 0. Hence, by setting q 2 = 



2 



-kvfvi 3 , thanks to (|21j) . we get: 



1 

„2 f3 

<n >=-^. (25) 
a 1 
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Equation (|7j) can be easily integrated obtaining < p >= pod^/a 3 . Finally, 
the independent equations ([5]) and ([9]) become, with /„ ^0 

* + 9/„(l-/.) " " + 3 ' (26) 

•• , / 2 (2/„-l) , a • 3g 2 / P l(l-/,) _ n 
^ + 2A(l-/,) + a /w 2^ -°" (27) 

Equations (|25p and (|26|) are the master equations of our model involving the 
scale factor a(t) at the scale of homogeneity and the void fraction f v {t). In 
the following subsection we obtain from (|26p - (|27p a function which allows 
to reduce the Buchert equations in a workable form. 

2.2 A useful function for the model 

First of all, equation ^ can be written in the following way 

a_ AttG qg 2 g A 

a--^-^ + 97^X) + 3- (28) 

By considering equations ([26]) . (p7|) and (f28j) we obtain 



6^ + 3^ - 2-i- - 6^—^— ^ = 3A. (29) 

a a 2 1-/, a (1 - /„) (1 - f v f 



(30) 



From equation (|20j) we get 

3(1 - /„)d 



3(1 - f v )a - f v a 
Moreover, we have 

i 

_ 8irGp a 3 _ g 2 f* -fi 
m ~ 3HW ' k ~ aW Q ~ 9f v (l-f v )Hf {6) 

Finally, by multiplying (f29|) by 2^GpQa^/[f v a — 3(1 — f v )a] and with the 
help of USD and l|2S |l -(|27 |l . we get 

dt{l-fj- (l-f v )a A - (62) 
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After integrating and using (|19p we have: 



(1 ~ 

(1 " fv) 



•m 



n F = e h i 




(33) 



where is a small integration constant (<< 1) which will be determined 
later. In the original FB model, in the limit A = 0, the function VLp is a 
constant. Its value is 1. In the FB model its presence is justified in the 
sense that, if we consider only matter in walls, then there must exist a con- 
stant critical density of reference (namely £If = 1)- This critical parameter 
denotes the density at the recombination era within walls. After the re- 
combination, the voids begin to expand and occupy the particle horizon. 
In our AFB model this picture is no longer valid. The function £lp begins 
with the value 1 at the recombination and it reaches asymptotically. We 
have denoted this function as a density parameter by recalling the original 
meaning of this function in the Wiltshire model. However, we are not able 
to indicate a clear physical interpretation for In our context, formula 
([33]) is very useful and allows to reduce, simplify and partially integrate the 
Buchert equations, as we see in the next subsection. 

2.3 Final form of the equations 

We now write down the Buchert equations in a workable form. Combining 
the equations ([26]) . ([27]) . ([30]) and ([33]) one obtains 



The main Buchert equations become (]26l) and ([34h together with the relation 
(]33p . We can further simplify these equations using manipulations similar 
to those found in [22]l23]. 

Namely, a first equation is obtained by multiplying equation ([26]) by /„ and 
subtracting equation ([34]) . A second equation is obtained by multiplying 
equation ([26]) by (1 — /„) and adding equation (]34"]) . As a result we get 



-(2/,-l) + 

n z 



2/^a _ <r 
3a a 




(34) 




(35) 




(36) 
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Thanks to (f3T]h equations ([35]) and ([36]) become: 



(37) 





A further manipulation can be done by multiplying equation (|37p by 



where u = f^jj^- Equations (|39p -(|40p represent the final form of our field 
equations. We are now ready to analyze the main features of our model. 

3 A study of the equations 

3.1 The case with A = 

In the limit A = we formally recover the solution in [22] for the Buchert 
equations (|39p and (|40p . but without clock effects. It is important to note 
that we do not use approximations. The absence of clock effects is only the 
outcome of our chosen gauge, i.e. the uniform time flow which is different 
from the one used by the original Wiltshire model. As a consequence, we 
formally obtain, in the limit A = 0, the same master equations of [22] 
[23] expressed in terms of the time t, but with a phenomenological lapse 
function set to unity (uniform time flow gauge). This happens because a 
different interpretation for the parameter /„, is present in our model due to 
the different gauge used. Hence, in this limit, the wall factor a w evolves 
exactly as an Einstein-de Sitter universe. Furthermore, as it is shown in 
[52~1 13"4] , the general A = solution possesses a tracking limit in which the 
factor a v evolves approximatively as a Milne universe at late times. For f v 
the tracking limit solution is given in [221 [35] : 



H Q x (l - f v ) 3 a and equation (J38D by a Q 1 H Q 1 /„ 3 a. We get: 




(39) 




(40) 



3f V0 H()t + (1 — / Vo )(2 + fv ) 



(41) 
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Formula (|4ip shows that in absence of a cosmological constant, the voids 
monotonically expand and asymptotically reach the value 1, i.e. they occupy 
the full volume horizon. The deceleration parameter q = 2(1 — f v ) 2 /(2 + f v ) 2 
is always positive and reaches asymptotically zero from positive values. 



3.2 The case with A ^ 

3.2.1 Late times behaviour of the solutions 

First of all, we consider the far future limit of our model. This behaviour can 
be exactly determined by analyzing the field equations. In this limit, both 
the wall density p w and the void density p v are negligible and asymptotically 
become zero. Hence, in this limit, both a w and driven only by the 

cosmological constant. As a result, for t — > oo we must have: 

a w (t) ~ a v {t) ~ eV^* -». a ~ eV^K (42) 

Expressions (|42|) states that asymptotically, since the universe is dominated 
by the cosmological constant, both walls and voids must evolve as a de Sitter 
spacetime. Thus, asymptotically we have 

H w (t -> oo) = H v (t -> oo) = H(t -> oo) = J~^. (43) 



As a result, by taking the asymptotic limit in (|26p . we obtain 

lim - /» - = (44) 
t->oc f v {\ - f v ) 

Hence f v (t — > oo) = independently on the behaviour of f v . Formula 
(|44p implies that the backreaction Q reaches asymptotically zero. At late 
times the Buchert equations for a w and a v decouple and the exact GR 
solutions dominated by the cosmological constant arise. Note that the same 
phenomenon appears in the original Wiltshire model, where at late times 
the solutions for a w and a v look like the exact GR solutions and Q vanishes. 
Moreover, from equation ()33|) . for the parameter tip, at early times we have 
dp ~ 1, while asymptotically £If — )■ 0. The exact expression for dp will be 
calculated later. 

Another interesting asymptotic limit is the one involving /„. Its value 
at early times, to be in agreement with WMAP data, is fixed to be a priori 
f„i « 1. From p51l or thanks to & and for f Voo we have 



^=1-^(1-4 (45) 
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where: 



?> 



lim -g- = k^. (46) 



Concerning k^, from equation (|39p we have 

1 fwi Q"ui® j w (1 



(47) 



a 3 

We multiply (|47|) by -4. By performing the limit t —> oo, we get 

koo = 7: r hm —2- 3. (48) 

Formula (|4"8"]h thanks to $32]), implies that fcoo 7^ 0. As a consequence, 
differently from the A = case, f Voo 7^ 1. For its determination, a numerical 
integration of the equations is necessary. In any case, some considerations 
can be made. In the limit A = 0, we must have f Voa = 1 and f Voa must be 
a function of f v i and of the adimensional quantity vAij. Since \/Aij << 1, 
on general grounds we expect that (1 — f Voo ) « 1. However, it is possible 
to calculate this value only by a direct numerical computation. We only 
mention this interesting difference with respect to the FB model. 

3.2.2 Explicit expressions for a w and £lp 

In this subsection, we integrate directly our field equations obtaining the 
exact expressions for a w and ilp which allow to obtain further asymptotic 
behaviours and confirm the ones above. After differentiating (I47j) and ap- 
plying again ([4"7]) we get the equation: 

2— + = A. (49) 

Equation ([4*9]) is nothing but the exact Einstein equation for a w . Hence, the 
wall factor a w solution of the Buchert equations is exactly the same of the 
ACDM model. This does not represent a surprise. In fact, in the model 
[221 [23] the wall metric evolves as the exact Einstein solution, while the void 
metric evolves in a different manner with respect to the decoupled exact 
Einstein solution. As a result we have 

2 




a w0 sinhs — — t . (50) 
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Moreover: 



cosh \\Mrt 



sinh | t 



Expressions (|50|) and (|51|) confirm the estimation 
more, we calculate £If- We obtain 



and 



(51) 



Further- 



si/ 



( 



cosh 



3A + 
4 % 



\ 



cosh | t 



(52) 



Equation (139]) can be easily integrated to give 



(1 - /„)3q = a 



^Om(l 
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OA 



^ sinh 3 



cosh 3 



3A, 



(53) 



Finally, for a w $ we have 



Qp / ^Om(l — £i 
f 3 - 



^OA 



cosh 3 



_3A f 

2 *"» 



(54) 



Equation (|40p cannot be integrated in a elementary way. For late times, this 
equation reduces to the exact Einstein equation, i.e. without backreaction. 
Note that we could gain advantage from this fact to obtain an approximate 
expression for f v (t). 

3.2.3 Early and late times behaviours for a v , f v and I w 

Equation (|40p can be useful to calculate the early times behaviour of our 
model. At early times, we have 



3AV 2 
—I ta, a v 



(55) 



where the behaviour for a v has been obtained by evaluating (|40p at t ~ tj. 
By putting again the above expression for a v in (|4"Uj) . we obtain 



a v o — oq 



9\ ^ H, 3 
4 



i (^Om^? 
f 3. 



(56) 
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After putting ([56]) in ([53]) . we get the relation f v i = ej. In this way, all the 
parameters of the model are specified. With the same technique, we can 



obtain the asymptotic behaviour of a v . We get a v = a Voo ev 3 where 

Concerning the asymptotic behaviour of /„, from (|39p at the leading term 
in t we have 

fv = f Voo ~ + o(l), (58) 

being B an integration constant. For B > 0, f v reaches its asymptotic value 
from below, conversely if B < 0. Since we expect (1 — f Voo ) « 1, perhaps 
B > 0. 

As an important consequence of ([31]) . we can study the asymptotic be- 
haviour of I w . In fact, from ([30]) we get 

(59) 



afv 
3(1- fv)a 



1 



Thus, both at early and asymptotic times I w = 1. This behaviour can be 
understood thanks to f)42[) . Furthermore, the behaviour given by (|59[) can 
potentially account for the broadly observed Hubble flow, without invoking 
clock effects. This concludes a preliminary study of our solutions. 

3.2.4 A formula for our model and its consequences 

We derive a very useful and intriguing formula containing all the cosmolog- 
ical parameters of our model. In fact, by combining the equations (|37|) and 
([38]) and evaluating them at a given present time to, we get 



1 — y ^0m(l — fv ) - ^ ^ + \ffn)\j ^Ofc + ^0A + ^0m ^1 + ^ ~ 

(60) 

where Qof, given by ([52]) . is also calculated at to- This equation is a gener- 
alization of a similar relation present in |22| . An important consequence of 
([60]) is provided by its domain 

f\)A > n 0m (-1 + ^J^) - n ok . (61) 
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The inequality (|6ip provides an absolute minimum for f^oA (at least for 
positive curvatures) in terms of the actual matter and curvature density 
parameters and of the age of the universe. Prom (|6ip and (|33p . we deduce 
that Ooa > fvi at the recombination era and reaches the asymptotic limit 
> (1 — fvoo) at late times i.e, thanks to (fTOj) . Q\ reaches asymptotically the 
value 1. 

As a first consideration, (JBTj) implies that, with Qom and f^ofc held fixed, a 
higher age of the universe requires a higher value for Ooa, i-e- more dark 
energy. Moreover, formula (I61D also says that, in order to mimic dark en- 
ergy, a large amount of negative curvature is needed. While a large negative 
curvature seems in disagreement with analysis of the CMB, this calculation 
was obtained within the ACDM paradigm and is not directly relevant here. 
It should be stressed that some studies |38[ [39] on the evolution of voids 
fraction in the ACDM context estimate an actual fraction of voids between 
(40 — 70)% of the total particle horizon volume. In these studies voids are 
identified as underdense regions. However, it should be noted that an un- 
derdense region may not have necessarily a negative curvature on average: 
a sufficient amount, for example, of dark matter can be enough to make the 
average curvature non-negative. 

In the following numerical examples, we assume that the fraction of voids 
present in the universe has negative curvature. We calculate the cosmo- 
logical parameters that allow to have (40 — 70)% of voids with hyperbolic 
curvature. We now give some crude numerical examples. As an example, we 
pose to ~ 14 Gyr with A ~ 10~ 35 /s 2 and so we have Qof — 1/(3.33) and we 
pose Qom = 0.3. Hence, from ([60]) . by setting (remember that to and flo m 
are held fixed) tlok = 0.1, Qq\ = 0.72, we have f VQ = 0.4 and Qoq = —0.12. 
By setting for example Qok = 0.2, f2pA = 0.72, QoQ = —0.22, we have 
f Vo = 0.6. Finally, by choosing Qq\ = 0.72, Uok = 0.28, QoQ = —0.3, we 
obtain f vo = 0.7. These crude estimations show that a large (observed) 
fraction of voids could be in agreement with a relatively small value for £lok- 
Furthermore, in the examples quoted above we have Qoq + £lok — 0. As 
a result, we could have an inhomogeneous model with a fraction of voids 
fvo — 0-4 in agreement with the one actually observed and with the same 
present day values for Ooa and S]o m of the ACDM model. This suggests that 
the backreaction Q, rather than acting as a dark energy, could act on void 
scales as an effective positive curvature, balancing the negative curvature of 
the voids themselves. By describing a cosmological model where 

fiom + OoA^l, ^os + ^ofc^O, (62) 
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and with a large fraction of voids, we can obtain a model mimicking the 
relation f^om + ^OA = 1 of the flat ACDM one with the further constraint 
Ooq + ^ofc = 0. 

As a final consideration for this section, note that inequality (|6ip can be 
seen also in terms of the parameter Oqq. Thanks to (|10|) we have 



and from (}6Tj) and ([63]) we deduce that the density parameter Qo m plays an 
important role. For an overview on the backreaction issue, see [40 ], l32l [33] . 

4 Conclusions 

Many attempts have been made in the literature in order to explain or rule 
out the cosmological constant in the present day cosmological models. In 
particular, within the Buchert formalism, the FB model [22, 23j has been 
recently presented to explain dark energy in terms of the gravitational en- 
ergy stored in voids. The amount of such a gravitational energy sufficient to 
rule out the presence of A is rather huge (~ 38% at present day). Within 
the formalism used for the FB model, it is possible to build a more realistic 
model of universe than the ones obtained with exact solutions as, for exam- 
ple, LTB metrics. In fact, to mimic dark energy, such exact solutions require 
a single huge void or hump of order of Gpc. Such a picture is in contra- 
diction with the known 'geography' of the universe, where a web structure 
seems to be more appropriate. An important issue is represented by the 
calculation of the effects of inhomogeneities on the value of A. It is very 
useful to build a non-perturbative model, within the cosmological constant 
paradigm, mimicking the present day web structure of the universe. Our 
main aim was to present a model allowing to study, in a non-perturbative 
framework, the effects of the observed inhomogeneities on the cosmological 
parameters and on the cosmic evolution. We have generalized the FB model 
by introducing a cosmological constant but without introducing ab initio 
clock effects. We have written the relevant equations in a simple way. In 
the limit A = the equations reduce to the form presented in [22] [23] . We 
partially integrate and study the asymptotic behaviour of the solutions. A 
first difference with respect to the FB model is that the cosmological con- 
stant inhibits the voids to fill asymptotically all the visible universe, namely 
fvoo < 1- Instead, an expected small fraction of space is asymptotically 
filled by walls. The amount of such a fraction cannot be deduced by an 




(63) 
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asymptotic study of the equations and can be obtained only by integrating 
the field equations (a numerical integration will be the subject of a separate 
paper). The expansion of the voids is different with respect to the FB model 
and could be in principle detectable. In fact, the role of the cosmological 
constant is to make the asymptotic behaviour of a(t),a w (t),a v (t) the same, 
contrary to the FB model, where the late time limit of a v (t)/a w (t) is diverg- 
ing. In particular, this phenomenon is visible in the behaviour of I w : this 
approaches the value 1 both at early and late times. The behaviour of I w is 
also interesting when analyzed in terms of the AFB model with clock effects. 
In fact, we can add clock effects of [23] in our model. In this case, the field 
equations for f v ,a(t),a w ,a v (t) remain unchanged, while the parameter I w 
becomes a phenomenological lapse function. As a consequence, also in the 
presence of A, clock effects are possible: time delay effects can coexist with 
a non-vanishing cosmological constant and we can also have cosmic clock 
effects which are interpreted in terms of quasi-local gravitational energy in 
presence of the mysterious A. A broadly uniform Hubble flow can thus be 
obtained with smaller clock effects than in [22 , 23J and a non-zero cosmolog- 
ical constant. This reasoning shows that so called Sandage-de Vaucouleurs 
paradox can also be solved with A / 0. 

As a further consideration, within our model we deduce a rather inter- 
esting formula relating all the cosmological parameters of the AFB model 
at a given present time. A remarkable consequence of this formula is the 
existence of a minimum actual value for Q\ (at least for positive curva- 
tures) expressed in terms of the age of the universe and the actual density 
parameters and O m . 

As a final consideration, note that our model predicts the way in which 
voids expand in presence of A. This behaviour could be object of a cosmo- 
logical test. In fact, methods to measure the expansion rate of voids are in 
preparation (see for example [41J). When these methods will be available, 
they will represent a way to distinguish between the inhomogeneous models 
present in the literature. 
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Appendix 



In our model, the observer is located within walls with zero curvature. The 
metric within the walls can be written as 

ds 2 w = -dt 2 + a 2 w [drjl + r]ldQ 2 ] (64) 

2 2 

j-t2 i (1 — /f) 3 ° r, 2 i 2 jr>21 

= -dt + 3 [(%, + r/ w dn \ 

f3. 

The hyperbolic voids metric is: 

ds 2 v = -dt 2 + al[dr]l+smh 2 (r] v )dn 2 ] (65) 

= -dt 2 + [drj 2 , + sinh 2 ^)^ 2 ] . 

f 3. 

J III 

It should be noticed that the scale factors a w , a v are not the ones obtained by 
solving the exact Einstein equations but instead they are obtained by solving 
the Buchert equations (j26|) and (j27|) with a non-vanishing backreaction. 
Concerning the metric at the scale of homogeneity with the scale factor a(t) 
given by (|16p we have: 

ds 2 = -dt 2 + a 2 dr] 2 + A(t, r])dn 2 , (66) 

where A(t,rj) is an area function satisfying 4w L Adr] = a 2 Vi(r]-}{), where 
rftl is the particle horizon radius. Note that also the metric (|66p is not an 
exact solution of Einstein equations, i.e. it is not a LTB metric, but rather it 
is obtained by solving the full Buchert equations (for more on the interpre- 
tation of (j66|) see [23 1 ) . The main astrophysical observations are by means of 
photons propagating along null geodesies. In the Buchert average scheme we 
do not average along the light cone. Nevertheless, we assume that one can 
describe the light cone with the averaged geometry (I66p . In practice, rather 
than averaging a bundle of null geodesies, we assume that the light cone is 
described by means of the fiducial averaged geometry (i66j) . Consequently, 
an observer located within walls, must relate his geometry given by (|64p in 
term of the Friedmann fiducial metric (I66p along the past null cone, i.e. the 
metric of the observer must be dressed (see [201 [21] ) by means of the Fried- 
mann bias given by (|66p . In contrast to many approaches to the Buchert 
equations which ignore the issue of the dressing of cosmological parameters 
and simply relate the volume average scale factor to the observed redshift 
[42 \ l4"3j 144] . we will adopt the matching procedure of Wiltshire |23j and 
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match the radial null sections of (|64p and (|66p . with the phenomenological 
lapse function set to unity in our case. We obtain: 

i 

d Vw = f ™ idV , ■ (67) 



(68) 



Thanks to (|67p . the wall geometry (|64p becomes: 
dsl, = -dt 2 + a 2 



drj 2 + (l-f v )lfj v 2 w dn 2 



The past null cone equation in (|68|) is given by: 

V= -• 69 

Jt a 

The distance-redshift relation di{z) in terms of the redshift z is given by: 

d L (z) = a (l + z)rj w , (70) 

^ = (1"A) 3 / — • 

■/t (l-/,)5a 

The functions f v , a are obviously obtained by solving the Buchert equations. 
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